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We discuss the classical and quantum properties of non-local modified Gauss-Bonnet gravity in 
de Sitter space, using its equivalent representation via string-inspired local scalar-Gauss-Bonnet 
gravity with a scalar potential. A classical, multiply de Sitter universe solution is found where 
one of the de Sitter phases corresponds to the primordial inflationary epoch, while the other de 
Sitter space solution — the one with the smallest Hubble rate — describes the late-time acceleration 
of our universe. A Chameleon scenario for the theory under investigation is developed, and it is 
successfully used to show that the theory complies with gravitational tests. An explicit expression 
for the one-loop effective action for this non-local modified Gauss-Bonnet gravity in the de Sitter 
space is obtained. It is argued that this effective action might be an important step towards the 
solution of the cosmological constant problem. 

PACS numbers: 



I. INTRODUCTION 

The discovery of the accelerated expansion of the late-time universe was the starting point of a wide 
spectrum of different theoretical constructions which aim is to provide a reasonable explanation of this accel- 
eration. The simplest qualitative possibility for such construction is to consider a gravitational modification 
of General Relativity (GR), as compared with the introduction of extra, exotic dark components of the 
energy in an ordinary GR scheme. There are a number of different candidates which qualify for gravita- 
tional alternatives of dark energy (a general review of those theories can be found in [II]). Recently, as a 
new key proposal for dark energy, non-local gravitational theories have been considered (JUL El- It has been 
demonstrated that some versions of these non-local gravities, which depend on the curvature and its (inverse) 
derivatives only, are definitely able to pass the Solar System tests @, Q • A scalar-tensor representation for 
such theories has been developed too [3J, making explicit the connection with local gravities. Furthermore, 
the possibility to construct a unified description of the early-time inflation epoch with the late-time accel- 
eration period becomes quite natural in such theories [|[ (for related, string-inspired non-local theories with 
scalars as dark-energy models, see also [5[ and the references therein). 
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In this paper we study classical and quantum aspects of the non-local Gauss-Bonnet gravities introduced 
by some of the present authors in in the de Sitter space. Their massive and scalar-potential versions 
are proposed, and their relation with local string-inspired scalar-Gauss-Bonnet gravity will be investigated. 
Using the above equivalence, the corresponding de Sitter solution with a constant scalar field will be explicitly 
constructed. The Chameleon scenario for this theory will be investigated too, and it will be shown that the 
theory under consideration passes the local tests (Newton's law, absence of instabilities). The one-loop 
effective action for the theory under discussion (again, using its equivalence with the local version) will be 
explicitly evaluated on the de Sitter space. The corresponding effective action is then explicitly obtained by 
using zeta-regularization, and it is used in the important discussion of the induced cosmological constant. 
Finally, some different version of massive non-local GB gravity, which may also be presented as a local 
multi-scalar-GB theory, is proposed. A classical, multiple de Sitter solution of this theory is found, where 
one of the de Sitter points can serve for the description of the inflationary epoch, while the other de Sitter 
solution, with much smaller value of the corresponding Hubble parameter, can be used for the description 
of the late-time acceleration period. 



II. NON-LOCAL MODIFIED GAUSS-BONNET GRAVITY AS STRING-INSPIRED 

SCALAR-GAUSS-BONNET THEORY 

We consider the non-local Gauss-Bonnet (GB) model introduced in Q, 

s = J*^(±-£gar*g) +Sm , (i) 

where S m is the matter action, R the scalar curvature, Q the Gauss-Bonnet invariant, and □ the 
d'Alembertian operator in the metric gij, with determinant g. Finally, k is related to the Newton con- 
stant G by k 2 = 8irG/c 3 . By introducing the scalar field 4>, one can rewrite the action (TTJ in a local form, 
namely 

S = J d A x^—g (2^-^ 9 ij dict>d j( t> + + S m . (2) 

2 -. 

In fact, one of the field equations of the latter action gives <\> — ~ ^-0~ L Q. By substituting this expression 
into ((2]), one obtains {T]). Note that the action J2J) corresponds to string- inspired scalar GB gravity, which 
was proposed as a dark energy model in Ref . [6|, [7J ■ 

For the sake of generality, we add a potential V(</>) to the Lagrangian density in ((2|), that is 

S = J (2^ ~ i ^di^djc), - V{ct>) + j>g\ + S m . (3) 

If we eliminate 4> by using the field equation \3<f> — V'(<p) — — Q, we obtain a rather untractable non-local 
theory (here and in the following V , V" mean derivatives with respect to the argument). There is however 
a simple case, that is 

V(<t>) = (4) 

where a mass term is present, which, typically, may be thought of as a non-perturbative string correction. 
In case of (|4]), by eliminating the scalar field (f>, we obtain 

s = J d i x^^-^g(D-m 2 y 1 g^j+s m . (5) 

Coming back to action ([3]), by assuming the metric to be the FRW one, with flat spatial section, and cf> to 
depend on the cosmological time only, the equations of motion — neglecting for simplicity the contribution 
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due to matter — read 
3H 2 



= — 



,.2 



+ ^~2<j> 2 + V(<P)-24<j)H 3 , 0=~^, + 3H^ -V'(<f>) + 24 (hH 2 + . (6) 



As usual, H — H(t) represents the Hubble parameter and the dot, as in H , </>, means derivative with respect 
to the cosmological time. If we further restrict to de Sitter space, that is, if we take H = Ho and <f> = <f>Q to 
be constant, we obtain 

V(4> ) = ^$, V'{M=2AHl. (7) 

K 

Thus, we see that there is always a de Sitter solution if the potential satisfies the condition {fj for a specific 
value of 4>- 

For example, if V is a mass term as in Eq. Q, then from ([7]) we get 

As a second — non trivial — example we consider the potential 

V = V e«* , (9) 
with positive constants Vq and go- The model has a de Sitter solution, where 

*-Hb&)- (10) 

In [J it has been shown that a classical de Sitter solution exists in the absence of the potential, too, but in 
such case the background field <f> is a time-dependent function. Extending this formulation one can construct 
de Sitter solution with time-dependent scalars also in the presence of the potential. However, such solutions 
will not be discussed here, due to fact that we are primarily interesting in the quantum properties of non-local 
modified GB gravity on the de Sitter background with constant scalars. 



III. CHAMELEON SCENARIO IN NON-LOCAL MODIFIED GAUSS-BONNET GRAVITY 

Many dark energy models generically include propagating scalar modes, which might render a large correc- 
tion to Newton's law, if the scalar field couples with usual matter. In order to avoid this problem, a scenario 
called the Chameleon mechanism has been proposed [8|, [9(. In this scenario, the mass of the scalar mode 
becomes large due to the coupling with matter or the scalar curvature in the Solar System and/or on the 
Earth. Since the range of the force mediated by the scalar field is given by the Compton length, if the mass 
is large enough-and therefore the Compton length is short enough-the correction to Newton's law becomes 
very small and it cannot be observed. The Chameleon mechanism has been used to obtain realistic models of 
F(i?)-gravity [HI, [H, EH . I n this section, we are going to consider the variant of the Chameleon mechanism 
which is generated by the coupling of the scalar field with the GB invariant in the action. A mechanism of 
this kind has been proposed in Q ■ As we have shown, the non-local GB theory can be rewritten as a model 
with a scalar field coupled with the Gauss-Bonnet invariant. Then, the Chameleon mechanism could work 
in a theory of this class, which we will investigate in the present section. 

For the action ([3]) the mass of the field is given by 

m = . 11 



Here <fi e is the background value of in a local region, like on the Earth or the entire Solar System. In the 
relevant region in which we are investigating the possible corrections to Newton's law, the curvature and the 
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background scalar field could be almost constant. Then the equation given by the variation of the action Q 
with respect to 4> takes the following form: 

V'(4>e) = Qe ■ (12) 

Here, Q e is the background value of Q. For example, on the Earth, we find that Q e ~ 10~ 71 eV 4 . 
For the model (g]) the mass ([TTj) is given by 

ml = . 13 

Since H ~ 10~ 33 eV and 1/k ~ 10 28 eV, by using Q we find 

-m-308 

ro! - — eV 2 . (14) 

v a 

In order that the Compton length could be 1 /im, that is, m<$, ~ 1 eV, we find that a ~ 10~ 308 , which is very 
small but, anyway, the Chameleon mechanism could perfectly work. 
For the model © Eq. (|T2j) has the form: 



and the mass (|TH)) is given by 



Then, by using (fTDj) . we find 



qV e"1" = G e , (15) 



a ^ 2 gT/ e^ 

mi = . (16) 



i rr 248 

m 2 -^ eV 2 . (17) 

v a 



Thus, when a ~ 10 248 the Compton length could be 1 fim. 

For a different example, we can now consider the following model, 

VW " iMTO ■ (18) 

Here Vo and /3 are positive constants, and a is a constant. As long as 4> > a, V(<fi) is a positive and smooth 
function of <j>. In this case, Eqs. © have the following form 

Vo _ 3# 2 V (30 2 - 2a0 o + /3 2 ) _ ^ ^ 



(0o-a)(^+/? 2 ) k 2 ' (0 o _ a ) 2 (02 +/32) 2 
By eliminating ffo in the two above equations (| 19[> . we obtain 



3</> 2 - 2a0 + /3 2 + -k 4 V = , (20) 



which can be solved as follows: 



q ± y/a 2 - 3/3 2 - 8« 2 Vfr 
<P = <P± = g ■ (21) 

Since <f>+ > as long as a < 0, we find c/>± > a. Then, V(<^>±) is surely positive and therefore there are 
two de Sitter solutions where the Hubble rate is given by 

Hl =h1 ^t L ± -am + (P)) ■ (22) 
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The smaller one could be identified with the present acceleratedly expanding universe and the larger one, 
with the inflation epoch in the early universe. We should note that, since there is no singularity in V(4>) for 
4>- < 4> < (f>+ , the two solutions are connected smoothly and, therefore, the transition from the inflationary 
era to the Dark Energy universe is possible, in principle. On Earth or at the Solar System scale, Eq. ([12")) 
acquires the following form: 

_ F o (30g-2a0 e+/ 3 2 ) = 
(0 e - a) 2 (<Pl + P 2 f 

which could be solved with respect to <j> e . The mass (fTB")) is given by 

_ n 2 Vg {I6a<t>l ~ 6 (a 2 + /? 2 ) - 2/3 2 (a 2 + /3 2 ) } 

We can conveniently choose the parameters so that rn 2 ^ becomes large enough in order not to give any 
measurable correction to Newton's law. 

We have thus shown that the Chameleon mechanism can actually work even in the new situation when we 
deal with a non-local GB theory which is equivalent to a local scalar-Einstein-GB theory, where the scalar 
field couples with the Gauss-Bonnet invariant. There, even though the non-local GB theory contains a scalar 
mode, this scalar mode does not provide any observable correction to the Newtonian law and, thus, a theory 
of this kind could emerge as a perfectly viable theory. Moreover, as we were able to see, the possibility to 
unify early-time inflation with late-time acceleration becomes quite natural in this context, what is an added 
bonus worth mentioning. Note also that the understanding of the equivalence principle in modified gravity 
[Io| may be somehow different from its standard formulation. 



IV. ONE-LOOP EFFECTIVE ACTION IN THE NON-LOCAL GAUSS-BONNET GRAVITY ON 

DE SITTER SPACE 

Here we discuss the one- loop quantization (for a review see Ref. [7J|] ) of the classical models we will be 
dealing with here, on a maximally symmetric space in the Euclidean approach. One-loop contributions can 
be important, especially during the inflationary phase. In any case, as it was shown in (lq . their analysis 
also provides an alternative method to study the stability with respect to non-homogeneous perturbations 
around de Sitter solutions in modified gravitational models, in agreement with [HI [17]. 

We start with the non-local GB-gravity related to the generalized model in ([3]). For the sake of simplicity, 
in the present section we will neglect the matter action S m , since it is irrelevant for our aims, and we shall 
use units in which the speed of light c = 1 and the Newton constant IQttG — 1. Our analysis is going to be 
very general: we will consider a non-minimal interacting term between the scalar and gravity of the form 
f(4>)G, with / being an arbitrary function. In this way, the model is described by the (Euclidean) action 

d 4 x^g [R + f(ct>)Q - ag ij d^djcf)- V {</>)] . (25) 

When / is a constant, this action reduces to Einstein's gravity minimally coupled to the scalar field. 

In accordance with the background field method, now we consider the small fluctuations of the fields 
around the de Sitter manifold (jjij, <j>o), of the kind 

R 

Rijrs = -Tg {{jirgjs - 9is9jr) , Rq = 4A = const , O = const . (26) 

This is a classical solution if the potential satisfies the conditions 
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which are the analog of (|7|) written in terms of the scalar curvature i?o- We see that, if V((f>) — ,then the 
Minkowski solution emerges. 

For the arbitrary solutions (gij,<j)) of the field equations, we set 



9ij = 9ij + h 



(28) 



and perform a Taylor expansion of the action around the de Sitter manifold, up to second order in the 
small perturbations (hij, if). Before we proceed with the expansion, it is convenient to write the action in 
order to take into account the fact that Q is a topological invariant. Then we observe that the non minimal 
interacting term between gravity and the scalar field, around the background solution, can be written in the 
form 



Gf(<t>) = Qf(<l>o)+G[f(<f>)-f(<f> )} 
= Gf(M + Go [/(<« - f(<t>o)} 



(G - Go) [/(<£) - /(&>)] 



(29) 



Here Go = Rq/® ^ s the value of the Gauss-Bonnet invariant evaluated on the de Sitter background. We note 
that the first term on the right hand side of the latter equation does not give contributions to the classical 
field equations and can be dropped out, while the second, proportional to Qq, modifies the scalar potential. 
Then, for our aim it is convenient to write the classical action (|25p in the final form 



d A x^ R + (/(0) - /(&>)) (G - g Q ) - agv d^d^ - Vtf) 



(30) 



where we have introduced the effective potential 

v(<t>) = F(0)-a o [/(0)-/(0 o )]. 



(31) 



One can check that action ([30)1 is equivalent to the original one ([2]) when /(</>) = <fi. 

We are now ready to perform the expansion. In the following we will use the compact notation Vq = V(cj>o), 
V Q ' = V'{<t>n), and so on. After a straightforward computation, along the same lines as for one-loop F(R)- 
gravity one obtains 



S[h] 



Ro-Vo-V<v+ ^-^)h + C 2 



(32) 



where £2 represents the quadratic contribution in the fluctuation fields (hij , tp) . Disregarding total deriva- 
tives, this reads 



/ V"\ 1 



f'Ro ( A + -i? ) + 2V ' 



i (2V0 - Ro 



-A- 



Ro 



& + h l - 



A ^| A rr 



(33) 



Here Vfc and A = § lJ VjVj represent the covariant derivative and the Laplace operator, respectively, in the 
unperturbed metric gij . We have also carried out the standard expansion of the tensor field hij in irreducible 
components [3, [l9[ , that is 



h^ = hj + Vi^j + Vj£i + VjVjcr+ -gij(h- Act) , 



(34) 



where a is the scalar component, while £j and hij are the vector and tensor components, respectively, with 
the properties 



Vif = , Vih) = , 



. 



(35) 



As is well known, invariance under diffcomorphisms renders the operator in the (h, a) sector non-invertible. 
One needs a gauge fixing term and a corresponding ghost compensating term. We consider the class of gauge 
conditions 

Xk = Vjhjk - —r- Vfc h , 



parameterized by the real parameter p. As gauge fixing, we choose the quite general term [1^ 

1 



£gf = § X* Gij x 3 , 



Gij =jgij +P A , 



(36) 



where the term proportional to 7 is the one normally used in Einstein's gravity. The corresponding ghost 
Lagrangian reads [lj| 



. Sx k 

£gh = B l Gik——C° 

5e J 



(37) 



where Ck and Bk are the ghost and anti-ghost vector fields, respectively, while 5 x h is the variation of the 
gauge condition due to an infinitesimal gauge transformation of the field. It reads 



5 h i:j = Vj£,' + VjEi 



j— = g i:j A + R tJ + — ^— ViVj • 



Neglecting total derivatives one has 



where we have set 



£ gh = B l (7 + /3 A Ha) C\ 



Iki=!h, { A ' t) ' ~^ V,V ' 



In terms of irreducible components, one finally obtains 

2 



C 



gf 



Ro 



3p 



Ro 



e A 1 + -^ ik + -^h A + v Aoa 



p 7 a 7 9 / . Rq , , 



2 



7? 



A QCT 



(38) 
(39) 
(40) 



(41) 



^gft = 7 



AnC 



Z>" I A , + ^ ) A 1 f 
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where ghost irreducible components are defined by 

C k — Ck+ VfeC , 
B k = B k +V k b, 



V fe C fe = 0, 
V k B k = 0. 



(43) 



In order to compute the one-loop contributions to the effective action we have to consider the path integral 
for the bilinear part 



C = C 2 + C 



gf 



C 



ah 



(44) 



of the total Lagrangian, and take into account the Jacobian due to the change of variables with respect to 
the original ones. In this way, the Euclidean one-loop partition function reads [Til . [T^ | 



= (det Gy)~ 1/2 J D[h i3 ]D[C k ]D[B k ] exp (- J d 4 xy/§£ 



(detGy)" 1/2 det Jf 1 det J 2 1/2 



£>[/i]D[ft ij -]£>[| J ']£>[<7]D[C fc ]£)[B fc ] J D[c]£>[6] exp 



d x \f]j C 



(45) 



where J\ and J 2 are the Jacobians coming from the change of variables in the ghost and tensor sectors, 
respectively. They read fbj ] 



Ji = A 0: 



Rc 



J 2 = -Ai— -Ao-^ A 



Ro 



Finally the determinant of the operator Gij , acting on vectors assumes, in our gauge, the form 



det d- 



const det I A i + ~ ) det I A + -fl + ~ 



7 











(46) 



(47) 



while it is trivial in the simplest case = 0. By A n (n — 0, 1, 2) we indicate the Laplacians acting on scalar, 
vector, and transverse tensor fields, respectively. 

Now, a straightforward computation, disregarding zero gravity modes and the multiplicative anomaly as 
well (|2y|)), leads to the expression of the one-loop contribution (7, j3, p) to the Euclidean partition 
function. It is a quite involved expression, which depends on the gauge parameters, and for this reason we 
will only write it explicitly in the Landau gauge corresponding to the choice 7 = 00, = 0, p = 1. On- 



shell, Z, 



,1) 

(7,/3,p) 



is independent of the gauge and it is compatible with a similar expression obtained in [19f | for 



Einstein's theory with a cosmological constant Ao = A = Rq/A, but in presence of a scalar field. In fact, we 
get 



shell 



det (-A 1 - A) 
det (-A 2 + § A) 



" 1/2 






det ^ 



V" 
1 



-1/2 



(48) 



The latter term is due to the scalar field, but it also depends on the coupling with the Gauss-Bonnet invariant. 
It has to be noted that the on-shell partition function is obtained by imposing conditions (|27| and Q = Rq/6. 
The latter condition in the expression of the gauge-dependent one-loop partition function is equivalent to 
dropping all terms proportional to /g. 

Now, we explicitly write the off-shell partition function in Landau's gauge. It reads 



Z, 



(i) 

(00,0,1) 



det (-A ! - A) 
det (-A 2 + §A- Vb) 



det -A 



-1 1/2 



det 



-An-^ 



Rq 

12 



q 1 det -A 



Rq 

12 



q 2 det -A 



Ro 
12 



<Z3 



-1/2 



(49) 
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The quantities qx, q 2 , (73, which depend in general on A, are the roots of the third-order algebraic equation 

q 3 + c 2 q 2 + Cl q + c = 0, (50) 



where 



Co 



4(/^) 2 ^ - 24/XVq' + 36(Vq') 2 - 36VoVy 



^ 2[10(/^) 2 ^ - 48/X^o + 36 Q E V-o + 54(U ') 2 - 9R V " - 18V V "} 



C2 = 



4[7/X + 9a(R Q + W Q ) - - 9U ")] 

R [6a+(f^Rl] 



The one-loop effective action can now be evaluated by making use of zeta-function regularization, and 
actually computing the zeta-functions C(s\L n ) related to the differential-elliptic Laplace-like operators L n 
PH . Using the same notations as in Ref . [l8[ , we have (see the Appendix) 



r (^,o,i) - 2 Qo UJ + 2 0l lTj~2 g2 lT"~ 



^Qo (7+91) -^Qo(7 + 92) -^Qo(7+93) , 





1 






-2 ( 


3o(| + ? 8 ) 



where 



,(a) = C' (0|i„//i 2 ) , L n = -A„-^(a-a n ), (53) 



C'(0|i„) = lim ^-C(*|in) =-logdeti n , (54) 

and a = 9/4, ai = 13/4, a 2 = 17/4. 

The effective equation for the induced cosmological constant can be obtained by varying the effective 
action with respect to A [l9[ . We would now like to study the role of the Gauss-Bonnet term. For the sake 
of simplicity, we choose a scalar potential satisfying the conditions V ' — 0, Vq" = 0. This means that, in 
the absence of the Q term, namely when /(</>) = 0, the potential is constant and the action becomes the 
Einsteinian one with a cosmological constant Ao = Vq/2, while, when /(</>) — 4>, there is a coupling of the 
form 4>(g - Go). 

We consider separately the two cases f(cj>) = and /(</>) = cj). In the first case we get qi = 6A0/A, while 
92 = <?3 = and so they do not give a contribution. Thus, we get 

rgU - § «, (f ) + \ 0. (f) - \ 0. (-f + %) - £ Ob (| + *) ■ (55) 

Here we have introduced the dimensionless variable y = Aq/A > 0, so the effective equation for the induced 
cosmological constant A can be obtained by taking the derivative of the effective action with respect to y, 
since 



dr di ar (1) 



dy dy dy 
where I is the classical action on the S A background, which reads 

967T 2 487T 2 An 48tt 2 



(00,0,1) =0; (56) 



1 = J d x^g(R -Vo) = j^— = -j— (2y-y ) . (57) 
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From (|A4|) and (|A5[) we get 



dQ n {a) 1 rM| f b k a k \, A , da ^ dF a (0) \ A 

\ i I I I ^ 1 _|_5 2 Q, ^g. 



<9y y \ £^ A:! / 3/i 2 y dy \ da J 3/x 2 y 




+ ai +761 + (a 2 +76 2 )a + VG(fc)a' £ 1 . (58) 



k=3 



For the first case, at lower order in y, from (|56[) we obtain the equation 

96tt 2 ,.. s 21.15 



A 



(1 - y) + — 54.60 + 41.32 y 

2/ 

+24 log -f - - 18 y log + 0(y 2 , y 2 log y) . (59) 
o/i y o/i y 



In the second case we also have vector and tensor contributions, as in the previous one but, in addition, we 
get three scalar contributions related with the roots of (|50"]) . Using (f52"|) and (j56"j) . at lowest order in y we 
obtain 

96tt 2 , x 20.83 / a \ 

— (1 - y) + 57.11 + 43.50 - 2.93 y 

A y V A o/ 

A o /, _ , nn a \ , A 



+251og-^-- 15-1.70-2 y log -^ + 0(y 2 ,y 2 logy). (60) 
o/x y \ A / ,5/1 y 

In Figs. |T]) and @ we have plotted the right hand sides of Eqs. (|5"9"j) and (f6"0"|) , respectively, as functions of 
y, for two different sets of values of the parameters. Notice that the zeros of these functions correspond, in 
each case, to a zero value of the induced cosmological constant. Aside from the cases here explicitly depicted, 
it can be easily seen that in our effective models the possibilities to get an induced cosmological constant 
which is exactly zero are reasonably high, since for values of the parameters lying in wide regions of the 
domain of expected values, either one or two roots of the equations exist, yielding the value of the induced 
cosmological constant exactly zero. As is well known [19(, there is a possible resolution of the the induced 
A-term problem coming from the contribution of higher-loop terms. This is explained in (l9l | precisely in 
the example of an effective action in pure Einstein gravity for a de Sitter background: one can get a very 
small effective A term irrespective of the tree level cosmological constant. But an even better possibility is 
to start from a zero tree level cosmological term, since quantum corrections will respect this property. We 
see that in our model there are good chances to realize this latter situation. One remark is in order. We 
got the one-loop effective action for non-local GB gravity using its classical equivalence with local scalar-GB 
gravity and working in terms of such local theory. It is quite well-known that such classical equivalence may 
be broken already at one-loop level. However, the equivalence is restored on-shell, i.e. using the one-loop 
corrected equations of motion. That is precisely the situation in which the induced cosmological constant 
has been here discussed. 

Now, we come back to the original action ([2]) with V{4>) — and f(<p) — <fi, which is equivalent to the 
non-local action ([T]). Also for this case there is a de Sitter solution gij,<fio, but cf>o is not a constant. This 
means that the term Gf{4>o) m (|29p cannot be dropped out and, thus, it gives a contribution to In 
principle, it is possible to take such contribution into account, but technically this is quite complicated, since 
it contains a lot of independent terms which mix scalar with vector and tensor components, so that the 
partition function is given in terms of the determinant of an involved 5x5 matrix of differential operators. 

When V((j>) = the original model has, however, a Minkowskian solution. This means that, using ([52"]) 
with V(4>) = 0, we can compute the cosmological constant induced by quantum fluctuations around the 
Minkowski background, and this may be interpreted as the spontaneous creation of a de Sitter universe 
starting from a flat one, which is a quite interesting feature. 
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FIG. 1: Plot of the right-hand-side of Eq. I)59[l vs y, for the particular values of the parameters Ao — fi 2 = 1 (left 
figure) and Ao = (J 2 = 10 (right figure). 





FIG. 2: Plot of the right-hand-side of Eq. (|60[1 vs y, for the particular values of the parameters Ao = (J 2 = a — 1 (left 
figure) and Ao = fi 2 = 10, a = 100 (right figure). 



V. OTHER CLASSICAL NON-LOCAL GB MODELS AND THEIR DE SITTER SOLUTIONS 

Another model giving rise to an interesting non-local action is the following 



S = I d xy/—g 



2k z 2a 



(61) 



where a is a dimensional constant and F(Q) is an adequate function of Q. By introducing three scalar fields, 
4>, £, and 77, one can rewrite the action (|6ip under the following form 



(62) 



s = J d^v^t^ + ^d^d^ + m^ + Hv-G^j ■ 

We may further add a potential V((f>) to the action 

s = J d^x^ ( A + ^d^d^ + 4>F{v) +Z(v-G)- V(4>) 



(63) 
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In the FRW universe, this action leads to the following equations 

= -\h 2 - -^0 2 + 24e& 3 + V(4>) , - 4 U + 3H<p) + F(tj) - V(4>) , 

= <j>F{rj) + £ , = r] - 24 (h 4 + HH 2 } . (64) 

However, after adding the potential V(4>), it is difficult to get the corresponding non-local action explicitly. 

When V(<f>) = 0, assuming that <fr = c^t, T) = ryo, £ = c^t, and H = Hq, with constant c^, 770, and i?o, 
the equations in (|64p reduce to the algebraic ones 

= -^H 2 - + 24c ? i/ 3 , = ^i/ o c + F(t7 ) , = c^fa,) + Cfi , = % - 24ff 4 . (65) 
We can solve Eqs. (|65|) with respect to c^, eg, and ?7o as follows: c<f, = ^ aH ' , = 



(24tf 4 ) F' (24ff 4 ), Va = 24iJ 4 , and we find 



= "ft - ( G o) 2 + ^G F (G„) F' (G ) . (66) 
8«r 18a 9a 



Here G = 24H$. 

For example, if we choose 

F{G)=hG 2 , (67) 

d^SJ) gives 

_1_ lk 2 / 2 N 
8k 2 + 18a "° I ' 



G o = 77^172 . ( 69 ) 



which has a trivial solution Go = corresponding to the flat background and 

9a 
44k*/, 

which corresponds to the de Sitter universe. 
As another example, one can choose 

F{Qf =-9o(Q+9i)(Q-92) ■ (70) 

Here go, g\, and gi are positive constants. In order that F (Q) can be real, we restrict the value of Q as 
-9i < G < 52- Then Eq. ([66]) yields 

= G 2 + -^G +.gi 52 , (71) 

which can be solved as 



1 / 9a / f 9a 2 



G0 = G ° ± "2| v -4^ ± Vl4^) +4/yi/y - I 

which are positive, provided a < 0, what we will assume in what follows. In order that —g\ < Gq± < 52, we 
find 
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Then, as long as Eq. (|73")1 is satisfied, there are two solutions describing a de Sitter universe. By choosing a 
more general F(Q), we may possibly find there can be several of these solutions corresponding to de Sitter 
universes, of which the solution with the largest value of Hq could perfectly correspond to the inflation epoch 
and the one with the smallest Hq to the accelerated expansion period of our present universe. Note that, in 
the model (|70[) . as far as Eq. ([73]) be satisfied, there is no singularity, say at Q — 0. Thus, a smooth transition 
could occur and there exists a true possibility to realize matter dominance before late-time acceleration. It 
is also possible to evaluate the one-loop effective action for the above theory, too, but the technical details 
of the calculation are quite involved in this case. 

VI. DISCUSSION AND CONCLUSIONS 

In summary, we have studied in this paper non-local GB gravity in its classically-equivalent local scalar-GB 
form, with a scalar potential. The classical de Sitter solution was obtained for the model with one scalar, as 
well as for the model with several scalars. In the last case, the multiple de Sitter solution encountered can 
be used for the unification of the early-time inflationary epoch with the late-time accelerating one (where 
an unstable de Sitter point should be considered for the inflationary epoch). The Chameleon scenario for 
non-local GB gravity has been presented and it has been shown that the theory discussed here satisfies the 
local gravitational tests (Newton's law) . 

The one-loop effective action for the theory under investigation has been calculated on the de Sitter 
background. Its explicit representation, in terms of a zeta-regularization scheme, gives us the possibility to 
determine, in all precision, the induced cosmological constant. It could well happen in this context that, as 
a result of the quantum corrections, such induced cosmological constant could eventually tend to zero, thus 
providing a natural solution to the elusive cosmological constant problem. 

As a first application of the results obtained in this work, it is quite interesting to observe that it becomes 
a rather immediate issue to investigate the stability of our model around the de Sitter solution. This is in 
fact necessary, in order to demonstrate that the de Sitter space is unstable, and most naturally provides a 
consistent, graceful exit from the inflationary era. However, that of taking into account infrared quantum 
gravity effects and the leading log approximation may eventually be a relevant issue, in order to construct 
the quantum gravity- induced inflationary universe (see e.g. [13|). To this aim, it is sufficient to require that 
the Laplace-like operators appearing in the on-shell, one-loop partition function (|48p be positive operators 
[l5| . As a result, we get 

Vq" = Vq" Go fS = Vo" ~ Rl fg>0. (74) 

For the linear case f((f) = <fi, this condition reduces to V " > 0, and this is certainly satisfied by the two 
potentials discussed in Sect. [TTJ 

The general, one-loop effective action on the de Sitter space we have found in this paper can be used as 
the basis for the search of the non-trivial RG fixed point in the functional, exact RG approach to string- 
inspired, scalar-GB gravity, in the same way as in scalar-Einstein/gauged SG theories [231] . with account to 
higher-derivative invariants [24j (in our case, the GB term). Finally, the one- loop calculation performed here 
is central to the subsequent investigation of the cosmological perturbation theory in non-local gravities. 

Acknowledgments 

We thank the referee of a previous version of this manuscript for remarks and criticisms that definitely 
contributed to its improvement. One of the authors (G. Cognola) acknowledges the support received from the 
European Science Foundation (ESF) for the activity entitled "New Trends and Applications of the Casimir 
Effect" (Exchange Grant 2262) for the period March-April 2009, during which the present work has been 
completed. He his also grateful to the Institute ICE/CSIC in Barcelona for the kind hospitality received 
during that period. This paper is, in part, also an outcome of the collaboration program INFN (Italy) - 



14 



DGICYT (Spain). It has been partly supported by MEC (Spain), projects FIS2006-02842 and PIE2007- 
50/023, by AGAUR (Generalitat de Catalunya), contract 2005SGR-00790, and by the Global COE Program 
of Nagoya University provided by the Japan Society for the Promotion of Science (G07). 



APPENDIX A: EVALUATION OF THE FUNCTIONAL DETERMINANTS 

We will here make use of zeta function regularization (see, for example, [2l|) for the evaluation of the 
functional determinants appearing in the one-loop effective action, Eq. (|4"5|) computed in the previous sec- 
tions. We shall first outline the standard technique based on binomial expansion, which relates the zeta- 
functions corresponding to the operators A, with eigenvalues A n > and A = j^(A — a), with eigenvalues 

A n = -§-(A„ — a), a being a real constant. With this choice, A„ and a are dimensionless. We assume to be 
dealing with a second-order differential operator on a D dimensional compact manifold. Then, by definition, 
for 3?s > D/2 one has 

COO = c(s|i) = £V. (Ai) 

n 

Us) = c(s\a) = J2 k s = (§) " £(a„ - «r s - (A2) 

where, as usual, zero eigenvalues have to be excluded from the sum. In order the binomial expansion in (|A2[) 
to be meaningful, we have to treat separately the several terms satisfying the condition |A„| > |a|. So, we 
write 



Us) 



Ro 
12 



k=0 



fcir(s) 



(A3) 



where we have set 



F a (s)= J2 CK~a)- s , F(s)= ]T A- s , (A4) 

A„<|a|; A„^a A„<|a| 

G(s)= V = C(s)-F(s), F(0)-F Q (0) =N , (A5) 

A„>|a| 



Nq being the number of zero- modes. It needs to be noted that (|A3[) is valid also in the presence of zero-modes 
or negative eigenvalues for the operator A. In many interesting cases, F a (s) and F(s) are vanishing and thus 
G(s) = C(s). 

As is well known, the zeta function may have simple poles on the real axis for s < D/2 but it will be 
always regular at the origin. Of course, the same analytic structure is also valid for the function G(s). One 
has 

r ^) = E B + (n _ m jD)/2 + ^). ( A6 ) 

71=0 y " 

J(s) being an analytic function and K n the heat-kernel coefficients which depend on geometrical invariants. 
In the physical applications we want to deal with, we have to work with the zeta function and its derivative 
at zero; thus, it is convenient to consider the Laurent expansion around s = of the functions 

T(s + k)C(s + k) = ^+a fe + 0( s ), (A7) 
s 

r(s + k)G(s + k) = —+a k + 0(s), (A8) 
s 
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bo 


= b k 


-F(0), 


ao = a +jF(0) , 


bk 


= h 


= K D -2k , 


a k = a k - T(k)F(k) , 


bk 


= b k 


= 0, 


G(k) = ((k) - F(k) , k> 



D 

1< k < — . 
~ ~ 2 



D 
~2 



(A9) 
(A10) 

(All) 



From previous considerations, one obtains 

<.<•)-(£ 



0<k<D/2 



b k a k (a k + -ib k )a k 



k\ 



+F a (s) + s 

k>D/2 



k\ 



a k G{k) 
k 



+ 0(s 2 ) 



and finally 



Ca(0) = F Q (0) 

C(0) - -Ca(0) 10. 



E 

0<fe<_D/2 

12 



b k a k 
k\ ' 

E 

0<fc<D/2 



(afc + 7bfc)a' i 
fc! 



+ ^(0)+ 2 

k>D/2 



a k G(k) 
k 



(A12) 

(A13) 
(A14) 



7 being the Euler-Mascheroni constant. If there are negative eigenvalues then F' a (0) has an imaginary part, 
which points out to an instability of the model. 

In the paper we have to deal with Laplace-like operators acting on scalar and constrained vector and 
tensor fields in a 4-dimensional dc Sitter space 50(4). In all such cases, the eigenvalues A„ and relative 
degeneracies g n can be written in the form 



A„ = (X - a) , g n = a (n + v) + c 3 (n + vf , A„ = (n + vf , 
where n = 0, 1, 2... and c\, c 2 , ^, a depend on the operator in question. In our case, we have 



L 



Li 



A ^° 

~ 12 9 



-A i q 

1 12 y 



3 

2 ' 



■ + <z, 



Cl 

1/ = 

Cl 

^ = 



Lo = -Ao - 



12 ' C 3 3 " 



C3 = l. 



T2 9 



1 12 ' 3 3 



(A15) 

(A16) 
(A17) 
(A18) 



where q are dimcnsionless parameters depending on the specific choice of f(R). 
We note that ((s) is related to the well known Hurwitz function (h (s, v) by 



oo oo 



2s-3 



n=0 



n=0 



and 



= Cl ( H (2s - 1, v) + c 3 ( H (2s - 3, v) 

G(s) = c 1 ( h (2s-1,v) + c 3 (h(2s-3,v)-F(s) 

= ciCh (2s -l, v + n) + c 3 ( H (2s - 3, v + h) , 



(A19) 



(A20) 
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n being the number of terms not satisfying the condition A„ > |a|. In order to proceed, we have to compute 
the quantities bk and a^, for k = 0, 1,2. To this aim, we note that the Hurwitz function has just a simple 
pole at 1, more precisely, 



1 



(A21) 



ip(s) being the logarithmic derivative of Euler's gamma function. After a straightforward computation, we 
get 



C(0)=ciCh(-1,^)+c 3 Ch(-3,^) 



bx = ~, & 2 = y, 



(A22) 



a = C'(0)-7C(0) 

= Cl [2& (-1, v) - 7 Ch (-1, v)\ + c 3 [2(' H (-3, !/) - 7 0f (-3, v)\ , 

71 



«i 



« 2 



"Cl 



ci (3, 2/) - c 3 



+ c 3 Cj?(-1,^) 

7-1 



Using (|A14|) . we obtain 

Q n (a) = C(0|i„/M 2 ) 



Fa(0) 



2 », 

■D 

fc=0 



fc! 



log 



i?0 

iV 



fc=0 



(at + jb k )a k 



fc=3 



a k G(k) 
k 



(A23) 
(A24) 

(A25) 



(A26) 



To conclude we would like to remark that in (| A4[) and (|A5|) it is not strictly necessary that A„ < a. In 
principle one can add up an arbitrary number of terms, with the only restriction that A„ ^ a, and all 
expressions we have derived here will be still valid. This means that n is allowed to be an arbitrary number. 
In this way, the convergence of the series in (|A26[) can be improved at will, which is a very nice feature of 
the procedure. 



[1] S. Nojiri and S. D. Od intsov, eConf C0602061, 06 (2006) [Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007)] 
[arXiv:hep-th/0601213| . 

[2] S. Deser and R. P. Woodard, Ph ys. Rev. Lett. 99, 111301 (200 7) [arXiv:0706.215l1 [astro-ph]]; C. Deffayet and 
R. P. Woodard. larXiv:0904.0961l [gr-qc]; N. A. Koshelev. arXiv:0809.4927 [gr-qc]. 

[3] S. Nojiri and S. D. Odintsov, Phys. Lett. B 659, 821 (2008) |arXiv:0708 .0924 [hep- th]]; S. Jhingan, S. Nojiri, 
S. D. Odintsov, M. Sami, I. Thongkool and S. Zerbini, Phys. Le tt. B 663, 424 (2008) |arXiv:0803.2613l [hep-th]] ; 
T. Koi visto, Phys. Rev. D 77, 123513 (2008) [arXiv:0803.3399l [gr-ac]] : T. S. Koivisto,~ Phys. Rev. D 78, 123505 
(2008) [arXiv:0807.3"778l [gr-qc]] . 

[4] S. Capozziello, E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Lett. B 671 (2009) 193 arXiv:0809.1535 [hep-th]]. 

[5] I. Y. Aref'eva, A. S. Koshelev and S. Y. Vern ov, Phys. Rev. D 72, 064017 (2005) [arXiv:astro-ph /0 507067] ; 
G. Calcagni, theory," JHEP 0605, 012 (2006 ) [arXiv:hep-th/0512259| ; I. Y. Aref'eva, L. V. Joukovskaya and 
S. Y. Vernov, JHEP 0707, 087 (2007) |arXiv:hep-th/0701184| ; G. Galcagni and G. Nardelli, larXiv:0904.4245l 
[hep-th]; G. Calcagni, M. Montobbio and G. Nardelli, Phys. Rev. D 76 (2007) 126001 arXiv:0705.3043] [hep-th]]. 

[6] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2005) [arXiv:hep -th/0504052 ; M. Sami, 
A. Toporensky, P. V. Tretjakov and S. Tsujikawa , Phys. Lett. B 619, 193 (2005); G. Calcagni, S. Tsujikawa and 
M. Sami, Class. Quant. Grav. 22, 39 77 (2005) [arXiv:hep-th/0505193] ; G. Calcagni, B. de Carlos and A. Do 
Felice, Nucl. Phys. B 752, 404 (2006) [arXiv:hep-th/060420i]; B. M. Leith and I. P. Neupane, JCAP 0705, 019 



17 



(2007) ; S. Nojiri, S. D. Odintsov and P. V. Tretyakov, Phys. Lett. B 651, 224 (2007) arXiv:0704.2520 [hep-th]]; 
S. Nojiri, S. D. Odintsov and M. Sami, Phys. Rev. D 7 4, 046004 (2006) [arXiv:hep-th/0605039j ; B. M. N. Carter 
and I. P. Neupane, JCAP 0606, 004 (2006) [arXiv:hep-th705 12262 ; T. Koivisto and D. F. Mota, Phys. Rev. 
D 75, 023518 (2007) [arXiv:hep-t h/0609155 ; G. Cogno la, E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, 
Phys. Rev. D 75, 086002 (2007) |arXiv:hep-th/0611198| ; M. R. Setare and E. N. Saridakis, Phys. Lett. B 670, 1 

(2008) ; A. K. Sanval. laxX~iv:0710.2450l [astro-phj ; B~C~Paul and S. Chose. lajXiv:0809.4131 1 [he p-th]; E. Elizalde, 
S. Jhingan, S. Nojiri, S. D. Odintsov, M . Sami and I. Tho ngkool, Eur. Phys. J. C 53, 447 (2008) |arXiv:0705.12TT1 
[hep-th]]; A. De Feli ce and T. Suyama , larXiv:0904.2092l [astro-ph.CO], 

[7] Y. Ito and S. Nojiri. larXiv:0904.0367l [hep-th]. 

[8] J. Khoury and A. Weltman, Phys. Rev. Lett. 93, 171104 (200 4) [arXiv:astro-ph/0309300] . 

[9] D. F. Mota and J. D. Barrow, Phys. Lett. B 581, 141 (2004) |arXiv:astrc-ph/0306047|~~^ 

[10] T. P. Sotiriou, V. Faraoni and S. Liberati, Int. J. Mod. Phys. D 17, 399 (2008)~ [arXiv:0707.2748"1 [gr-qc]]. 

[11] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007). 

[12] S. Nojiri and S. D. Odintsov, Phys. Lett. B 657, 238 (2007) [arXTv:0707.1941l [hep-th]]; S. Nojiri and 
S. D. Odintsov, Phys. Rev. D 77, 026007 (2008) [arXiv:0710. 17381 [hep- th]]; G. Cog nola, E. Elizalde, S. No- 
jiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, ibid. 77, 046009 (2008); [arXiv:0712.40T7l [hep-th]]. 

[13] S. A. Appleby and R. A. Battye, Phys. Lett. B 654, 7 (2007); L. Pogosian and A. Silvestri, Phys. Rev. D 
77, 023503 (2008); S. Nojiri and S. D. Odintsov, Phys. Lett. B 652, 343 (2007) [arXiv: 0706. 13781 [hep-th]]; 
S. Capozziello an d S. Tsujikawa, Phys. Rev. D 77, 107501 (2008); S. Tsujikawa, Phys.~ Rev. D 77, 023507 (2008) 
|arXiv:0709.T39"T1 [astro-ph] ] . 

[14] I. L. Buchbinder, S. D. Odintsov and I. L. Shapiro, Effective action in quantum gravity, IOP Publishing, Bristol, 
1992. 

[15] G. Cognola and S. Zerbini, J. Phys. A 39, 6245 (2006) [arXiv:hep-th/0511233| . 
[16] V. Faraoni, Phys. Rev. D 72, 124005 (2005) [arXiv:gr-qc/0511094| . 

[17] G. Cognola, M. Gastaldi and S. Zerbini, Int. J. Theor. Phys. 47, 898 (2008) [arXiv:gr-qc/0701138]. 

[18] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, JCAP 0502, 010 (2005) [arXiv:hep-th/056 l096 . 

[19] E. S. Fradkin and A. A. Tseytlin, Nucl. Phys. B 234, (1984) 472. 

[20] E. Elizalde, L. Vanzo and S. Zerbini, Commun. Math. Phys. 194 (1998) 613 [arXiv:hep-th /9701060] . 

[21] E. Elizalde, S. D. Odintsov, A. Romeo, A. A. Bytsenko and S. Zerbini Zeta regularization techniques with 

applications, World Scientific, 1994; E. Elizalde, Ten physical applications of spectral zeta functions, Lecture 

Notes in Physics, Springer- Verlag, Berlin, 1995; A. A. Bytsenko, G. Cognola, L. Vanzo and S. Zerbini, Phys. 

Rept. 266 (1996) 1 arXi v:hep-th/950506l] ; A. A. Bytsenko, G. Cognola, E. Elizalde, V. Moretti and S. Zerbini, 

Analytic aspects of quantum fields, World Scientific, Singapore, 2003. 
[22] N. C. Tsamis and R. P. Woodard, Class. Quant. Grav. 26, 105006 (2009) [arXiv:0807. 5006 [gr-qc]]. 
[23] D. Benedetti, P. F. Machado and F. Saueressig, larXiv:0902.4630l [hep-th]; L. N. Granda and S. D. Odintsov, 

Phys. Lett. B 409, 206 (1997) [arXiv:hep-th /9706062| . 

[24] A. Codello, R. Percacci and C. Rahmede, Annals Phys. 324, 414 (2009) |arXiv:0805.290 9 [hep-th]]. 



